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Abstract — In this paper the geometric process is used 
for the analysis of accelerated life testing for Marshall- 
Olkin Extended Exponential (MOEE) distribution using 
Type II censored data. Assuming that the lifetimes under 
increasing stress levels form a geometric process, The 
parameters are estimated by using the maximum 
likelihood method and the original parameters instead of 
the developing inference for the parameters of the log 
linear link function are used. The asymptotic interval 
estimates of the parameters of the distribution using 
Fisher information matrix are also obtained. The 
simulation study is conducted to illustrate the statistical 
properties of the parameters and the confidence intervals. 

Keywords — Maximum Likelihood Estimation; Type II 
Censoring; Survival Function; Fisher Information 
Matrix; Asymptotic Confidence Interval; Simulation 
Study. 

1. Introduction 

Models and methods of accelerated life-testing are useful 
when technical systems under test tend to have long 
lifetimes. Under normal operating conditions, as systems 
usually last long, the corresponding life-tests become too 
time-consuming and expensive. In these cases, accelerated 
life tests (ALT) can be applied to reduce the experimental 
time and hence the cost. There are basically three types of 
accelerated life tests: constant-stress test, step-stress test 
and progressive stress. 

In constant-stress test each experimental unit put at only 
one of the stress levels and in step-stress the level of stress 
is increased step by step until all items have failed or the 
test stops for other reasons. Progressive-stress loading is 
quite like the step stress testing with the difference that the 
stress level increases continuously 

Failure data obtained from ALT can be divided into two 
categories: complete (all failure data are available) or 
censored (some of failure data are missing). Due to 
different types of censoring, censored data can be divided 
into time-censored (or type I censored) data and failure- 
censored (or type II censored) data. Time censored (or type 
I censored) data is usually obtained when censoring time is 
fixed, and then the number of failures in that fixed time is a 
random variable. Failure censored (or type II censored) data 
is obtained when the test is terminated after a specified 
number of failures, and then time to obtain that fixed 
number of failures is a random variable. For more details 
about ALTs one can consult Bagdonavicius and Nikulin 
[i], Meeker and Escobar [ii], Nelson [iii, iv], Mann and 
Singpurwalla [v]. 

Constant stress ALT with different types of data and test 
planning has been studied by many authors. Watkins and 
John [vi] considers constant stress accelerated life tests 
based on Weibull distributions with constant shape and a 



log-linear link between scale and the stress factor which is 
terminated by a Type-II censoring regime at one of the stress 
levels. Ding et al. [vii] dealt with Weibull distribution to obtain 
accelerated life test sampling plans under type I progressive 
interval censoring with random removals. Pan et al. [viii] 
proposed a bivariate constant stress accelerated degradation 
test model by assuming that the copula parameter is a function 
of the stress level that can be described by a logistic function. 
Chen et al. [ix] discuss the optimal design of multiple stresses 
constant accelerated life test plan on non-rectangle test region. 
Fan and Yu [x] discuss the reliability analysis of the constant 
stress accelerated life tests when a parameter in the generalized 
gamma lifetime distribution is linear in the stress level. Islam 
and Ahmad [xi] and Ahmad and Islam [xii] discuss the optimal 
constant stress accelerated life test designs under periodic 
inspection and Type-I censoring. Geometric process (GP) is 
first used by Lam [xiii] in the study of repair replacement 
problem. Kamal [xiv] estimates Weibull parameters in 
accelerated life testing using geometric process with type-II 
censored data. 

This article is to focus on the maximum likelihood method for 
estimating the acceleration factor and the parameters of 
Marshall-Olkin Extended distribution. This work was 
conducted for CSALT with type II censored scheme. The 
confidence intervals for parameters are also obtained by using 
the asymptotic properties of normal distribution. In the last, 
the statistical properties of estimates and confidence intervals 
are examined through a simulation study. 

2. Model Description 
2.1 The geometric Process 

A GP is a stochastic process {X n ,n = 1,2,... } such that 

{X n ~ 1 X n ,n = 1,2,...} forms a renewal process where X>0 is 
real valued and called the ratio of the GP. It is easy to show 
that if {X n ,n = l,2,...} is a GP and the probability density 

function of X l is f{x) with mean \i and variance a 2 then 

the probability density function of X n will be A." _1 /(A." _1 x) 

with mean \ilX n ~ x and variance a 2 /A. 2( "~ 1) . It is clear to see 
that a GP is stochastically increasing if 0 < X < 1 and 
stochastically decreasing if X > 1 . Therefore, GP is a natural 
approach to analyze the data from a series of events with trend. 

2.2 The Marshall-Olkin Extended Exponential distribution 

(MOEE) 

Marshall and Olkin [xv] proposed a new method for adding a 
parameter to a family of distributions. Suppose we have a 
given distributionwith survival function 

(SF) F(x),— oo < x < oo then the Marshall-Olkin extended 
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distribution is defined by the SF 

— (xF(x^) 

G(x) = = — -co<x<oo,a>0,a = l-a (1) 

1 - aF(x) 

If the survival function of exponential distribution 
is F{x) = e~ Bx , x, 9 > 0 , and put it in equation (1), we 
obtain the SF 

— 1-a _ 

G(x) = — = , -oo<x<oo, a>0, a = l- a (2) 

e x - a 

The distribution with the survival function (2) is called the 
MOEE with parameters a and 9 . The probability density 
function (pdf) and the cumulative distribution function 
(cdf) and the hazard rate of the Marshall-Olkin extended 
exponential with the survival function (2), respectively are 
given by 

a9e 9x 

s(x;a, 8) = — q ^, -oo<x<oo,a>0,a = l-a (3) 

[e Bx -af 

e Qx -l 

G(x;a,9) = ^- (4) 

r( ~ x ) = ^7 (5) 

[« -a] 

when a = 1, the pdf, cdf, SF and hazard rate reduce to 
those of the exponential distribution. 

Assumptions 

• Let there are s increasing stress levels and under 
each stress level n items are inspected. Let the 
random variable x k ^ denote the lifetime of the 

tested product, where k, k = 1,2, s index for 

stress level and i,i = 1,2 n is the index for 

test item under each stress. In the Type II 
censoring scheme, the test is terminated at each 
stress level when r failures are observed. Under 

the k' h stress level, only the r smallest 
observations of the sample will be collected, 
whose failure time can be written 

as x km < x k{2 ) ^ ^ x k ^ r) . 

• The product life follows a Marshall-Olkin 
Extended Exponential distribution given by (1) at 
any stress. 

• Let the sequence of the random 

variables X 0 ,X l ,X 2 , ,X n , denote the 

lifetimes under each stress level, where X 0 

denotes item's lifetime under the design stress at 
which items will operate ordinarily. We assume 

{x k ,k = 1,2,3, s} is a geometric process with 

ratio X > 0 . 

• At any constant stress level S k , the mean life 
1/9^. is a log linear function of stress, that is, 
logllA^ ) = a + bS k , where a and b are unknown 
parameters that depend on the nature of the 
product and the test method. When k = 0 , the 
above equation depicts the relation of the mean 
life and the designed stress level. 
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2.3 G.P assuption in ALT 

From the above assumption, it can be easily shown that 



log 



-b(S k+] -S k ) = bAS 



Above relation can be written a 



-bAS 



It tells that when the increased stress levels form an arithmetic 
sequence with a constant difference AS , the mean life under 
each stress level forms a geometric sequence with the 
ratio e- hAS . 

Let e~ bAS = X , it is clear from the above relation 
fl A- = ^k-i = ^~Qk-2 = = X k Q 

In case of MOEE distribution, the pdf of a test item at the 
k' h stress level is: 

g x (x\a, Q,X) 



a9,A r aM k _ ie ™*-i x , t aQe e}}x 



r 6 k x — n2 r w k-\ x 

[e * -a] [e k 1 



a 



(e 6Xx -a) 2 



This implies that 

g Xk {x)=X k gxo (X k x) 

Therefore the pdf of a test item at the k stress level is 



8x, {x)=X k g(X k x) = X k 



uQe 



6A. x 



(e Qxkx -af 



(6) 



Consequently, the CDF of the test item at the k' stress level 
is: 



G x (x\a,Q,X)-- 



er.v 



1 



The probability that an item censored at time x k( r , is: 
1-a 



Gy (x)- 

k{r) V I 



8k K x kr — 

e h - a 



(7) 



(8) 



For those r k items that fail at time t , their order statistics can 
be denoted by X k ^ with PDF: 



1,2,. 



(9) 



3. Maximum Likelihood Estimation 

Let the test terminates at r' h failure time x^^and (n - r)units 
are still surviving at each stress level. The likelihood function 



of the k' h stress can be expressed as above 



8x 



■■8x k ( X k(r))[ G X k { X Hr)t ' 



k (n — r) 
(10) 

We derive the ML estimates of a, 6, and A from the 
likelihood function given by (6). The substitution of (6) and (8) 
in (10) gives 
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(n - r). 



it 

;=i 



a9e 



0/. . 



er. 



-l + a 



\e h +a-l. 



.0<x k[l) < <x k(r) (ll) 

It follows that the likelihood function of observed data in a 
total S stress levels is: 

f \ 



n 



n) 



(n - r). 



if 



OA.' j,, 



-l + a 



-a-l 



0<^ (l )< <x k{r) ;l<k< S .{\2) 

The log -likelihood function corresponding (12) can be 
rewritten as 



+ S ' (h-j 
4=1 

MLE's of a , 9 and X are obtained by solving the equations 

— = 0, — = 0 and — = 0 , where 

3a 59 dX 



ln a — lnl Xjir + a - 1 









1 + aj 



>ZZ 



l 



dl ns 



da a ^— ' I ex' x u , . — 

4=1 /=! e 41 + a-l 4=1 



c e».**, +a _j 



(13) 

dl ra.v^vr,* 



9 ZX'" H ' "ZZ 

4=1 !'=1 



er. 



. ; . = i \e~ -l + a 
{n-r)k x kr e J - 



k=i \e eXx *'+a-l 



(14) 



dl 



s r s M t-i ex**., 



— i — r J9fcA, 



z 

it=l 



,ea** t , 



+ a-l 



(15) 



The MLEs of a, 9, ant/sexist but do not have a closed 
form. The Newton Iterative Method is applied to 
obtain d, 9, and X . 

4. Fisher Information Matrix 

The Fisher information matrix F is obtained by taking the 
negative second partial derivatives of the log -likelihood 
function and can be written as 



d 2 l 


d 2 l 


d 2 l 


da 2 


5a59 


dadX 


d 2 l 


_dH_ 


d 2 l 


595a 


59 2 


dQdX 


5 2 / 


d 2 i 


dl 


dXda 


DXdQ 


dX 2 



Elements of Fisher Information matrix are 



8 I ns 

a 2 ~ 2 
da a 



d A l _ rs 

ae 2 e 2 



-zz 



k=l 1=1 



-z 



-1+af **0*^+a-lf 



2 ZZ 



/ ,\,2t 2 ex*x t , 

(a-lJA. x, e i: 



-l + a 



(n- rfa-ljk x kr e k < 
e Ar + a-l 



5A, 2 



k 



2 IX 



ZZ>ZZ^-^ 

=1 1=1 4=1 i-1 

9*(ife-l)X, M j: H e m *- ,: « +(a-l)|e*(*:-l)A.''" 2 x (:i e e ^-' t '' + e 2 4 2 A M ~ 2 x 2 .e ia '' I * i 



2 2 If 

a 2 / e 2 / , 



3aee aeea 



£Z 

4=1 i=i 



9 9 s r 

a 2 / a 2 / = £ 



dadX dkda 



=zz 

4=1 i=l 



e 6X - x " -l + a 



e b -l+a 



z 

4=1 



-z 

4=1 



{n-r)X k x kr e eXx *r 



k >- -l + a 



e *' -l + a 



-2 



zz 

4=1 i=l 



4=1 i = l 



e 9X ^ -l + a 



4=1 



tr - 1 + ( 



5. Asymptotic confidence interval estimates 

Now the variance and covariance matrix of parameters can be 
written as 
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d 2 l 

da 2 
d 2 l 

dQda 
d 2 l 

dXda 



d 2 l 

dadQ 
d 2 l 

ae 2 

d 2 l 

dim 



d 2 l 



dadX 
d 2 l 

dk 2 



AVar(a) 
ACovjed) 
ACo\Akd) 



AVai{d) 
ACoii^d) 



A CovldX, J 
ACoJfii) 
AVar\k) 



The 100(l - p)% asymptotic confidence interval for a, 9 
and X are then given respectively as 



d + Z p ^AVar(a) 



e+z 



2 



and 



X + Z p^AW^k) 



6. Simulation Study 

The performance of the estimates can be evaluated some 
measures of accuracy which are the mean squared error 
(MSE) and the coverage rate of asymptotic confidence 
intervals for different sample sizes and stress levels. The 
main purpose of simulation is to obtain MLEs of 
estimator and judge their performance for this model. 
To perform the simulation study first a random sample 

x ki ,k=l,2, ?, / = 1,2. ...r is generated from MOEE 

distribution which is censored at r = 30,35 . 
The value of the parameters and numbers of the stress 
levels are chosen to be 0 = 0.2 , a = 0.9 , X = 4 and 
5 = 2,4. 

For different sample sizes n= 40, 60, 80,100 & 200 
ML estimates, the mean squared error (MSE), lower 
and upper asymptotic confidence interval limits and 
the coverage rate of the 95% confidence interval of 
parameters are obtained. The results obtained in the 
above simulation study are summarized in Table 1, 2, 
3 &4. 

Table I 

Simulations Results Based on Censored Data from 
GP MOEE with 9=0.2, a =0.9 and X =4 for s=2 & 
r=30 



n 


e 

d 


MSE (e) 
MSE (d) 
MSE (i) 


95% 
Confidence 
Interval 


95% 
Confidence 
Interval 












i 




LCL 


UCL 


Coverage 




1.1134 


0.3213 


1.1577 


1.1456 


0.94342 


40 


0.0524 


0.0464 


0.0504 


0.0634 


0.94569 




1.1426 


0.7464 


1.1345 


1.1789 


0.94433 




1.1078 


0.2345 


1.1096 


1.2535 


0.95896 


60 


0.0456 


0.0345 


0.0367 


0.0486 


0.96789 




1.0565 


0.6234 


1.0455 


1.0765 


0.95662 




1.1024 


0.2245 


1.1018 


1.1065 


0.95644 


80 


0.0367 


0.0312 


0.0254 


0.0465 


0.96547 




1.0534 


0.4424 


1.0341 


1.0655 


0.95778 




1.0234 


0.2145 


1.0214 


1.0356 


0.96654 


100 


0.0354 


0.0234 


0.0343 


0.0455 


0.97986 




1.0134 


0.4134 


1.0123 


1.0256 


0.98766 




1.0154 


0.1575 


1.0145 


1.0167 


0.96775 


200 


0.0165 


0.0187 


0.0157 


0.0169 


0.97997 




1.0121 


0.3567 


1.0112 


1.0322 


0.98979 
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Table 2 

Simulations Results Based on Censored Data from GP 
MOEE with 9=0.2, a =0.9 and X =4 for s=2 & r=35 



n 


§ 

d 
X 


MSE (e) 
MSE (d) 
MSE [X ) 


95% 
Confidence 
Interval 


95% 
Confidence 
Interval 
Coverage 


LCL 


UCL 


40 


1 9811 
1 .Zo 1 1 

0.0635 

1.5734 


yj.j 1 j 1 
0.0774 
1.3837 


1.1179 
0.0573 
1.2772 


1.3645 
0.0874 
1.6677 


0 ens'* 1 
0.92984 
0.93784 


60 


1.2534 
0.0624 
1.5646 


0.5717 
0.0676 
1.2847 


1.1763 
0.0563 
1.3636 


1.3771 
0.0776 
1.7731 


0.93736 
0.95341 
0.93857 


80 


1.1837 
0.0731 
1.5663 


0.3733 
0.0635 
1.2841 


1.1254 
0.0582 
1.4883 


1.2764 
0.0767 
1.6554 


0.95736 
0.93836 
0.93726 


100 


1.1532 
0.0701 
1.4525 


0.2736 
0.0573 
1.1836 


1.1235 
0.0563 
1.3746 


1.3746 
0.0781 
1.6831 


0.95837 
0.94736 
0.95836 


200 


1.0328 
0.0572 
1.4281 


0.4786 
0.0463 
1.1577 


1.0273 
0.0471 
1.3745 


1.2774 
0.0589 
1.7663 


0.95732 
0.96722 
0.96721 



Table 3 

Simulations Results Based on Censored Data from GP 
MOEE with 9=0.2, a =0.9 and X =4 for s=4 & r=30 







MSE 


95% 






e 


(8) 

MSE 


Confidence 
Interval 


95% 


n 








Confidence 


d 


(a) 






Interval 




X 


MSE 


LCL 


UCL 


Coverage 




1.1255 


0.2654 


1.1087 


1.2787 


0.93567 


40 


0.0543 


0.0656 


0.0503 


0.0697 


0.93786 




1.6445 


0.2997 


1.6043 


1.7675 


0.94886 




1.1144 


0.2575 


1.1129 


1.1209 


0.93864 


60 


0.0524 


0.0624 


0.0487 


0.0579 


0.94785 




1.5770 


0.2769 


1.3458 


1.7807 


0.94324 




1.0865 


0.2522 


1.0346 


1.2543 


0.95464 


80 


0.0475 


0.0586 


0.0267 


0.0586 


0.94456 




1.4764 


0.2345 


1.3467 


1.5578 


0.94342 




1.0531 


0.2235 


1.0154 


1.0757 


0.95888 


100 


0.0374 


0.0478 


0.0235 


0.0389 


0.95890 




1.3664 


0.2257 


1.1686 


1.4689 


0.96868 




1.0134 


0.1755 


1.0127 


1.0174 


0.96454 


200 


0.0167 


0.0165 


0.0158 


0.0198 


0.97675 




1.2565 


0.2142 


1.0458 


1.2896 


0.97865 
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Table 4 

Simulations Results Based on Censored Data from GP 
MOEE with 9=0.2, a =0.9 and X =4 for s=4 & r=35 



n 


9 

i 


IVlOL* 

(e) 

MSE 

(a) 

MSE 

w 


95% 
Confidence 
Interval 


95% 
Confidence 
Interval 
Coverage 


LCL 


UCL 


40 


1.1134 
0.0524 
1.1426 


0.3213 
0.0464 
0.7464 


1.1577 
0.0504 
1.1354 


1.1456 
0.0634 
1.1789 


0.94342 
0.94569 
0.94433 


60 


1.1078 
0.0456 
1.0565 


0.2345 
0.0345 
0.6234 


1.1096 
0.0367 
1.0455 


1.2535 
0.0486 
1.0765 


0.95896 
0.96789 
0.95662 


80 


1.1024 
0.0367 
1.0534 


0.2245 
0.0132 
0.4424 


1.1018 
0.0254 
1.0341 


1.1065 
0.0465 
1.0655 


0.95644 
0.96547 
0.95778 


100 


1.0234 
0.0354 
1.0134 


0.2145 
0.0235 
0.4143 


1.0214 
0.0243 
1.0123 


1.0356 
0.0455 
1.0256 


0.96654 
0.97986 
0.98766 


200 


1.0154 
0.0165 
1.0121 


0.1575 
0.0187 
0.3567 


1.0145 
0.0157 
1.0112 


1.0167 
0.0169 
1.0322 


0.96775 
0.97997 
0.98979 



7. Conclusions 

From the results in above tables, it is easy to find that 
estimates of the parameter perform well. For fixed 9, a 
and A., the MSEs decreases as n increases and the 
probability coverage is also satisfactory. This indicates that 
the ML estimates provide asymptotically normally 
distributed and consistent estimator for the parameters. For 
the fixed sample sizes, as number of failures r gets larger 
the MSEs of the estimators decrease. This is very usual 
because more failures increase the efficiency of the 
estimators. 

From above discussion and results it is concluded that the 
present model work well under type-II censored data for 
MOEE distribution and would be a good choice to be 
considered in the field of ALTs in future. For the further 
research in this direction one can choose some other 
lifetime distribution with different types of censoring 
schemes. 
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